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Abstract 



We study the rest-frame instant form of a new formulation of relativistic 
perfect fluids in terms of new generalized Eulerian configuration coordinates. 
After the separation of the relativistic center of mass from the relative vari- 
ables on the Wigner hyper-planes, we define orientational and shape variables 
for the fluid, viewed as a relativistic extended deformable body, by introducing 
dynamical body frames. Finally we define Dixon's multipoles for the fluid. 
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I. INTRODUCTION 



Both theoretical and numerical investigations concerning relativistic hydrodynamics are 
becoming very important for astrophysics, cosmology and also for heavy-ions collisions. 
Therefore it is important to develop analytical methods able to describe the various aspects 
of the theory of relativistic perfect fluids, in particular their properties as isolated extended 
relativistic systems. 

Usually relativistic fluids are described by assigning i) a unit four- velocity field U^{z) 
\U^Un — 1]; ii) the local thermodynamical functions internal energy p{z), particle density 
n{z), local pressure p{z), local temperature T{z), entropy per particle s{z)); iii) an equation 
of state. Here z^^ are the coordinates of a fixed point in Minkowski space-time and some- 
times this description is named Eulerian. However some clarification is needed about this 
terminology at the relativistic level, where a notion of simultaneity has to be introduced 
to give sense to the relativistic equations of motion, which very often are deduced from an 
action principle. 

In the standard approach to non-relativistic fiuid dynamics, not based on variational 
principles, two distinct points of view are usually used ^: the Lagrangian (or material) point 
of view and the Eulerian (or local) point of view. The equations of motion, i) the continuity 
equation or mass conservation; ii) Euler-Newton equations or balance of linear momentum; 
iii) the conservation of energy; look different in the two points of view. 

A) In the Lagrangian point of view the fluid is described by the flux lines x(t, Xo) with 
x{0,Xo) — Xo, defined as the integral lines of the 3- velocity field u{t,Xo) — ^2^^. Each 
integral line is labeled with its initial coordinate Xo- The coordinates the Lagrangian 

(or comoving) coordinates of the Lagrangian point of view. The fiux lines have the role to 
describe the mechanical aspect of the fiow of the fiuid. If we think that associated with each 
fiux line there is a material particle defined by an infinitesimal volume of fiuid around the 
point Xo, then the flux line x{t, Xo) is also the trajectory followed by the material particle, so 
that also the name material point of view is used. For a flxed value of Xo, x{t, Xo) specifles the 
path of the mass element which was at at t = 0; for a flxed value of t, x{t, Xo) determines 
the transformation of the region initially occupied by the whole mass of the fluid. For every 



See for instance Section 3.2 of Ref. [1]. 



3 



local thermodynamical function evaluated on the flux lines Q{t, x{t, Xo)) its expression in the 
Lagrangian point of view is Q{t,Xo) — Q{t,x{t,Xo))- Therefore the equations of motion in 
the Lagrangian point of view are written using only the ordinary time derivative, ^ Q{t, Xo)- 

B) Instead in the Eulerian (or local) point of view the fluid is described by taking as 
Eulerian coordinates a set of coordinates x referring to a fixed location in space and not 
to a moving mass element of fluid and by using the 3- velocity field v{t,x), which is the 
velocity of the fluid particle that happens to be at the location x at time t. If at time t we 
make the identification x — x{t,Xo), namely the fixed coordinate is seen as the coordinate 
of the flux line through that point, we recover the 3- velocity field of the other point of view: 
v{t,x)\^=^(t,xo) = u{t,Xo). Now a local thermodynamical function is described by the local 
function Q{t,x) with G(t,x)\£=£(^t,xo) = G{t,Xo), so that this description of the fluid is also 
named the local point of view. The equations of motion in the Eulerian point of view involve 
the so-called total (or material) derivative -^Q^tjx) = {§[ + u{t,x) ■ Q{t, x), since we 
have I g{t, Xo) = g g{t, x{t, fj) = (| + u{t, Xo) ■ -§g) Q{t, x{t, f^)) = g{t, f)) \s=x{t,xo)- 

To extend these descriptions to the relativistic level in the Minkowski space-time with 
coordinates z^^ without introducing an explicit breaking of covariance like the one implied 
by the decomposition z'^ = {z° = ct;z), it is convenient to work in the context of Dirac's 
parametrized Minkowski theories [2,3] on arbitrary (simultaneity and Cauchy) space- like 
hyper-surfaces, leaves of the foliation associated to a 3+1 splitting of Minkowski space-time. 
If T is the scalar parameter (mathematical time) which labels the leaves of the foliation, 
a are curvilinear coordinates on them (with respect to an arbitrary centroid x'^{r) — z'^{r, 0) 
chosen as origin) and z^{T,a) are the embeddings of the hyper-surfaces E,- in Minkowski 
space-time, then every local thermodynamical function has an equal time re-formulation: 
Q{z{T,a)) = Q{T,a). By using these adapted coordinates we have the following natural 
relativistic generalization of the two points of view. 

A) In the Lagrangian (or comoving) point of view the fluid is described by the flux lines 
(^{zo,f), defined as the integral lines of the four- velocity field with initial condition z^ and 
parametrized by their proper time f 

^C'i^cf) = u^mzo,f)), cn^o,o) = (1.1) 

The points used for the initial conditions must belong to a space-like Cauchy hyper- 
surface in Minkowski space-time. If we use the embeddings z^{t, a) of the hyper-surfaces E^ 
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of a foliation, the flux world-lines are described by functions C,'^{t^Bo) — C'^{zo,t^{t, ^o)) — 
z''{r,f{r,ao)) with f(0, ct<,) = 0, S(0, ct<,) = ('"{'^,^0) = zi'{ro,ao) = z^. Since r is not 
the proper time of any flux hne, on we have ^^^^^ / ^{^^^ = C/'^(r, Bo) for the flux 
line through ao at r = 0. 

In these adapted coordinates the coordinates Bo on the Cauchy surface S^-^ are the 
Lagrangian (or comoving) coordinates of the Lagrangian point of view ^ replacing the non- 
relativistic Xo, while the functions S(r, (Jo), describing the flux lines, replace the non- 
relativistic x{t,Xo) and u{t, Bq) — ^^^q^°^ is the 3- velocity field replacing u{t,Xo). The local 
thermodynamical functions evaluated along the flux lines Q{C^{t, Bo)) — G{z^{t, E(t, Bo))) — 
^(r, E(t, CTo)) have the expression ^(r, CTq) = ^)\ff=^{T,ao) ^ Si'^^'^i'T: ^o))- The equa- 
tions of motion involve only the ordinary time derivative ^ ^(r. Bo). To each flux line there 
is associated a material particle at Bo on the Cauchy surface St-=o- 

B) In the Eulerian point of view the fluid is described by taking the fixed 3-coordinates B 
(replacing the non-relativistic x) as Eulerian coordinates and r as the time in the coordina- 
tization dictated by the embedding. The 3- velocity field is v{t, B) and it can be shown that 
its expression in terms of the 4- velocity field U^{t,B), written in adapted coordinates, is 

— * 

v{t,B) — U{t,B)/U'^{t,B) [see Eq.(3.16)]. When the Eulerian coordinate is identified with 
the fiux line (''(r. Bo) = z^'{T, E(t, Bo)) we get v{t, o')U=s{r,ff„) = ^^fe^- The local thermo- 
dynamical functions are described by local functions Q{t, B) with Q{t, E(t, Bo)) = Q{t, Bg). 
Their equations of motion in the Eulerian point of view involve a total r-derivative, 

^ Oir, B) = i^^+ v{t, B) ■ ^) ^(r, a), since ^ ^(r, a,) = + • ^) ^(r, S(r, Bo)) = 

{^Q{r, ^))\ff=EiT,aoy 

Therefore, while in the Lagrangian (or material) point of view we follow the evolution 
of a thermodynamic function evaluated in a material particle (namely along the fiux line 
physically determined by the average particle motion), in the Eulerian (or local) point of 
view we follow the evolution of the same function evaluated in a given space-time point 
coinciding with the associated material particle only on the Cauchy surface. 

Till now only in the non-relativistic framework of the Euler-Newton equations there has 
been a study of the transformation from Eulerian to Lagrangian coordinates [4]. 

The next problem is how to derive the equations of motion of the Lagrangian and Eulerian 
point of views as the Euler-Lagrange (EL) equations of an action principle. 

An important development in relativistic hydrodynamics has been given by Brown [5], 
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who built a general framework encompassing all the known variational principles for rela- 
tivistic perfect fluids and clarifying the inter-connections among very different Lagrangian 
approaches. 

In Brown's paper the fluid is described by means of a set of scalar 3-dimensional configu- 
ration variables (5'(z), i = 1, 2, 3, of the space-time coordinates z^, interpreted as Lagrangian 
(or comoving) configuration coordinates for the fluid labeling the fluid flow lines. 

In Ref. [6] one of the action principles of Ref. [5] has been re-formulated in the contest of 
Dirac's parametrized Minkowski theories [2,3] on arbitrary (simultaneity and Cauchy) space- 
like hyper-surfaces, leaves of the foliation associated to a 3-1-1 splitting of Minkowski space- 
time for arbitrary equations of state of the type p{z) — p{n{z), s{z)). Now the Lagrangian 
(or comoving) coordinates of the fluid are q;*(t, a) — a^{z{T, a)) ^. For each value of r, we 
can invert a' = Q;*(t, a) to a — a{T, a') and use the q;*'s as a special coordinate system on 
inside the fluid support Vair) C E,-: 2;''(t, ^(t, a*)) = z^{T,a^). This approach is reviewed 
in Section II. 

Since it can be shown [sec Eqs.(2.16) and(3.2)] that we have (^^{t, ff)\ff=j2{T a ) ~ 
a*(0,(To) = <^o(^o)) namely that the coordinates a*(r, a) are constant along the flux lines, 
this explains why these coordinates are a possible set of Lagrangian (comoving) coordinates 
for the fluid in alternative to the ^o's . Therefore every action principle studied in Ref. 
[5] generates Euler-Lagrange equations [hke Eq.(2.19)] which describe what happens at the 
flxed location a when r changes, i.e. these EL equations generate the equations of motion 
of the thermodynamical functions Q{T,a) in the Eulerian (or local) point of view. This 
is a consequence of the necessity that the conflguration variables of an action principle be 
T-dependent, even when, like in this case, they are used to simulate the flxed comoving 
coordinates of the fluid. 

In this paper we will show that it is possible to deflne a variational approach whose conflg- 

— * 

uration variables are the adapted 3-coordinates E(r, (To) of the flux lines, which evolve from 

^The fluid is supposed to have compact support la(r) C S,-, whose boundary 5V^(r) is dynami- 
cally determined as the 2-dimensional surface in each 11^ where the pressure vanishes, p{z{t, a)) = 
for a G dVair). In the case of N disjoint fluid sectors, we can use the same description with 
Va{T) = Uj Vai{T) till when the compact supports Vaiir) do not overlap. 



6 



the Lagrangian coordinates Bo — S(0, on E^=0) instead of the Lagrangian (or comov- 
ing) variables a(T, ct) used in Ref. [5]. The configuration 3-coordinates E(r, CTo) are strictly 
speaking neither Lagrangian nor Eulerian coordinates. However they can be considered as 
generalized Eulerian configuration coordinates, because through the position a = S(r, cTq) 
(connecting the two points of view) they allow to obtain the Eulerian description of what 
happens in the fixed point ffo when r changes. Indeed in this case the resulting EL equations 
(3.17) [replacing the Eulerian ones (2.19)] will correspond to the equations of motion of the 
thermodynamical functions Q (r, ap) in the Lagrangian (or material) point of view. It will 
be shown (see footnote 7) how we can transform these EL equations into the equations of 
motion in the Eulerian point of view for Q (r, a) . This will allow to obtain a Hamiltonian 
formulation with Eulerian coordinates using Poisson brackets instead of the Lie-Poisson 
brackets of Ref. [7] . 

We shall study the Hamiltonian first class constraints of the fiuid on arbitrary space-like 
hyper-surfaces and, then, the restriction to the rest frame foliation (Wigner hyper- planes) , 
in order to obtain its rest-frame instant form, already used for particles and fields in Refs. 
[8-12]. Then, following the methods of Refs. [13-15], we shall study the problem of the 
center-of-mass and relative variables, the separation of relative variables in orientational 
and vibrational ones by means of the introduction of dynamical body frames and Dixon's 
multipoles [16,17] of the fiuid. 

In this way we get a complete control on the relativistic kinematics of perfect fiuids 
considered as extended deformable objects. The next step will be to couple the perfect fiuid 
to metric and tetrad canonical gravity, whose rest-frame instant form has been developed 
in Refs. [18-20], both to develop a scheme of Hamiltonian numerical gravity in accord with 
constraint theory and to study the linearized theory in a completely fixed 3-orthogonal 
Hamiltonian gauge following the scheme of Ref. [21]. Another future development [22] 
will be to study the non-relativistic limit of this approach and, then, after the addition of 
Newton gravitational potential, to recover the ellipsoidal equilibrium configurations [23] in 
this kinematics in the case of incompressible fiuids. 

While in Section II we review the description of perfect fiuids with Lagrangian coordi- 
nates, in Section 111 we introduce the new formulation with Eulerian coordinates. In Section 
IV we derive the Hamiltonian formulation associated with the action written in the previous 
Section and we get the usual constraints of Dirac's parametrized Minkowski theories and 
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their Dirac Poisson algebra [2] . In Section V the rest frame instant form of the dynamic 
is constructed. This form of the dynamic is such that we can discuss the problem of the 
separation of the relative variables from the center of mass-like variables and the analogous 
problem for the rotational and shape variables in the same way as it has already been done 
for relativistic particles. This is done in Sections VI and VII. In Section VIII we discuss the 
various either exact or approximate forms in which the invariant mass of the fluid, i.e. the 
Hamiltonian in the rest-frame instant form, may be presented as a function of the orienta- 
tional and shape variables. Dixon's multipoles for the fluid are defined in Section IX. In the 
final Section there are some concluding remarks. 

Appendix A reviews notations on space-like hyper-surfaces. In Appendix B there is a list 
of the equations of state for which we can obtain a closed form of the fluid invariant mass. 
Appendix C contains remarks on Poisson brackets. Appendix D describes the Gartenhaus- 
Schwartz transformation. Finally in Appendix E there are some solutions for the kernels 
associated with the relative variables. 
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II. LAGRANGIAN FORMULATION OF THE DYNAMICS OF RELATIVISTIC 
PERFECT FLUIDS IN PARAMETRIZED MINKOWSKI SPACE-TIME WITH 
LAGRANGIAN COMOVING CONFIGURATION COORDINATES. 



In this Section we review some of the results of Refs. [5,6]. One of the many action 
principles for a Lagrangian description of relativistic perfect fluid dynamics described in Ref. 
[5] has been re-formulated in Ref. [6] in the context of the parametrized Minkowski theories 
[8,3]. As said in the Introduction the starting point of these theories [2] is the foliation of the 
Minkowski space-time by a family of space-like hyper-surfaces defined by the embedding 
z^{t, a) {R xT, ^ M^) ^. The fields z^{t, a) define a coordinates transformation z^ i— > a"^ 
on the pseudo-Riemannian manifold M^, a field of cotetrads {/i are the flat indices) 



and the induced metric 



See Appendix A for other properties of space-like hyper-surfaces. 

In parametrized Minkowski theories the Lagrangian of every isolated system is written 
as a functional of the Lagrangian coordinates «*(r, cr) of the system, adapted to the folia- 
tion, and of the embedding z^{t, a) interpreted as the Lagrangian coordinates describing the 
hyper-surface in this enlarged configuration space. This functional is determined by consid- 
ering the Lagrangian of the system coupled to an external gravitational field and replacing 
the 4- metric g^ui^) with the induced metric (2.2) in the adapted coordinates. 

As said in the Introduction in the Eulerian point of view, the relativistic perfect fluid is 

characterized by the 4- velocity field U^{z) {U^U^ = 1) on Minkowski space-time and by 
a set of local thermodynamical functions. After the foliation of with the hyper-surfaces 
Et, the 4- velocity field has the adapted covariant components 

^Wc use the notation cr"^ = (t, cr') where A = (r, f), f = 1,2,3, for these coordinates adapted 
to the foliation. The notation r will be reserved in Section V for the 3-vectors on the Wigner 
hyper-planes 




(2.2) 
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(2.3) 



With the adopted parametrization of Minkowski space-time the local thermodynamical 
functions can be seen as functions of (r, a) by means of the replacement — z^{t, a) ^ 
. In particular let us consider the numerical density of particles n{T,a). Together with 
the 4- velocity field U^{t, a) it defines the numerical density current n{T, a) U^{t, a). The 
conservation of the total particle number is the following constraint on this current (" ;" and 
" ," denote the covariant and ordinary derivative, respectively) 

1 d 



n(r, a) U^{t, a) 



,a) da^ 



^Jg{r, a) n{T, a) U\t, a) 



= 0. 



(2.4) 



We can also to consider the energy density function p{T,a), the local pressure p{T,a), 
the local temperature T{t, a), the entropy per particles s(r, a) and the chemical potential 

Mr,g) = '''^- (2.5) 

n{T, a) 

The first principle of thermodynamics is given by the following differential relation 

dp{T, a) = ij,{t, a) dn{T, a) + n(r, a) T{t, a) ds{T, a). (2.6) 



The equation of state is given interpreting p as a function of n and s 

p(T,a) ^ p(n(T,a),s(T,a)), 



(2.7) 



and we can obtain the other thermodynamical quantities as functions of n and s, in particular 

1 d 

T{t, a) = — — ^ Tr(^('^> s(t, a)), 



n{T, a) ds 

p{t, a) = n{T, a) ^{n{r, a), s(t, a)) - p{n{T, a), s(t, a)). 
Finally, we have to add the entropy conservation 



(2.8) 



(2.9) 



This constraint tells us that we don't have loss of entropy out the flux tube defined by the 
fluid's flux hues. 



^Prom now on we shall denote with /(r, a) the functions /(r, a) = f{z{T, a)). 
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Due to the constraints (2. 4), (2. 9), the four-velocity field U^{T,a) and the independent 
thermodynamic functions, n{T,a), s{T,a) are a set of redundant variables for the fiuid de- 
scription . In Ref. [5] it is shown that the constraints (2. 4), (2. 9) may be enforced introducing 
some Lagrangian (or comoving) variables a^r, a) for describing the fluid. It is useful to in- 
troduce the flelds J'^(t, a) 

^giT,a)niT,a)u\T,a) = j\T,a). (2.10) 

These flelds are dependent on the Lagrangian coordinates q;*(t, a), i = 1, 2, 3, according to 
the deflnition 

J^{r,a) = -det I{T,a), 

r(r, a) = det I(r, a) [r\r, a)]', MLi^I , (2.11) 

where the 3x3 matrix 7(r, a) is 

r / ^-..i ( daHr, a)\ , , 

[/(T,^)r. = (^^Hj- (2-12) 

From the deflnition (2.11) it follows that (cr'^ = r) 

^Ar,a) = 0. (2.13) 

This equation is equivalent to the constraint (2.4). In particular, it follows from Eq. (2.13) 
that the total number of particles on the hyper-surfaces E^- may be deflned as {Va{T) is the 
fluid's volume on E,-) 

Af= f d^aJUT,a), (2.14) 

and it is conserved by the evolution in the time parameter r. In this parametrization the 
entropy per particle is a function of the q;'(t, a) alone 

s = s{a\T,a)). (2.15) 

By construction it satisfy the entropy constraint (2.9) [see also the following equation (2.16)]. 

The Lagrangian coordinates tt*(r, c?) can be interpreted as resulting from a coordinate 
transformation o"'' i-^ a* on the hyper-surface E^. In particular a*(0, cr) define a coordinate 
transformation cr'' i— > a\0,a) on the hyper-surface E^=o; if Ki(0) is the fiuid total volume 
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in this hyper-surface, every point in the total volume on the hyper-surface E^, V^(t), is in 
a one to one correspondence with a point in Va{0) by means of the flux lines. Due to the 
definitions (2.10) and (2.11) it follows that the fields q;*(t, ct) are constant along the flux 
lines, since by construction we have 

u\T,a)^a'{T,a)^0. (2.16) 

Then the fields Q;*(r, a) can be interpreted also as labels assigned to the flux lines; the field 
q;'(t, a) tells us that the flux line labeled with a* goes through the point z'^{r, a) e M^. 

With the previous definitions and observations, the action defined in Section 5 of Ref. 
[5] has been rewritten in Ref. [6] in the form 

S^-J drd^a^g{T,a)p{n[a],s[a]), (2.17) 

where s is given by Eq.(2.15) and from Eq.(2.10) it follows that 



n(r, a) = ^ — ■ (2-18) 

The stationarity of the action with respect to variations of the a's gives the fluid equations 
of motion ^ 

= 2 V^E,F]i^, a) U^\t, a) + T(r, a) 5,^(r, a), (2.19) 

where 

VA{T,a)^f,{T,a)U^{T,a), (2.20) 

is the Taub vector, if //(r, a) is the chemical potential (2.5). In the previous relation we have 
used the notation 

V v ' - YMZ^ZA 

^[E,F] — ^[E;F] — 2 ■ l^-^-l-J 

The variation of the action with respect to the metric variations Sgj^j^ defines the stress- 
energy tensor 



^They are equations of motion in the Eulerian (or local) point of view. 
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5S 



p{t, a) +p(t, a) U\t, a) U^{t, a) - p{t, a) ^^^(r, a). 



(2.22) 



In Ref. [5] it is showed that the equations of motion (2.19) are equivalent to the stress- 
energy tensor conservation law 



To see this, we have to observe that Eq.(2.23) is equivalent to 



BC 



(2.23) 



{9AB - UaUb) T^"- -n = 0. 



(2.24) 



The first of these equations is equivalent to the entropy constraint (2.9), which is then 
satisfied. 



The second equation is equivalent to the Euler equations 



(2.25) 



(p + p) Ua.,b + {-^1 + U") P,B = 



(2.26) 



The stationarity of the action (2.17) with respect to the variations of the -2'^(r, a)'s gives 
us the equations of motion of the embeddings z^{t, a). Since the action depends on z^{t, a) 
only through the induced metric, they are 



5S 



d 



SS 



■ Ze>[T,a 



d 



S9Ab{'^:^) ^ 
\f^)T^\r,a)zUT,a) 



= Vy(r,a)^J(r, a) Qcoir, a) [r^^(r, a 



. - 0. 



(2.27) 



Due to the stress-energy conservation (2.23), following from the fluid equations of mo- 
tion (2.19), these equations arc always satisfied without any restriction on the embeddings 
z^{t, a), which remain arbitrary. In other words, the equations of motion (2.19) and (2.27) 
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are not independent. This is the Lagrangian manifestation that the parametrized Minkowski 
theories are singular theories. In these theories the z^{T,a) are gauge variables and in the 
Hamiltonian formulation their conjugate momenta are defined by first class Dirac constraints 
[3,8] as it will shown in Section IV . 
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III. LAGRANGIAN FORMULATION OF THE DYNAMICS OF RELATIVISTIC 
PERFECT FLUIDS IN PARAMETRIZED MINKOWSKI SPACE-TIME WITH 
EULERIAN CONFIGURATION COORDINATES. 



In this Section we introduce a different parametrization of the action (2.17) using a 
new set of configuration coordinates. The geometrical interpretation of the old Lagrangian 
(comoving) coordinates a*(r, a) as labels for the fiux lines suggests that it is possible to 
parametrize the action (2.17) with new adapted 3-coordinates E(r, CTq), which describe the 
flux fines as the integral curves of the 4-velocity field starting from the Lagrangian coor- 

— * 

dinates (To = S(0,CTo). As explained in the Introduction they can be named generalized 

— * 

Eulerian configuration coordinates. The derived variables C'*(''') ^o) — z^i^, S(t, CTo)) define 
the four-dimensional flux line going through ao si t — 



(It 



e(r,f?o) 



= C/'^(C(T,a„)), (3.1) 



with the initial condition ^^(0, CTq) = z^{0, Bo), namely E(0, Bo) — Bg. 
Moreover, due to Eq.(2.16), we get consistently 

a\T, S(r, Bo)) = a*(0, Bo) = <(ao), (3.2) 

i.e. the functions al{Bo) (a possible set of Lagrangian coordinates replacing the Bo^s) are 
constant along the fiux line through each point Bo- Then the inverse function theorem 
impfies E(r, Bo) — F{t, a(0, Bo)) and Bo — E(0, Bo) — -F(0, a(0, Bo)). Therefore at r = the 
generalized Eulerian coordinates of the fiuid are just the Lagrangian 3-coordinates Bo of the 
points where the fiux lines intersect the hyper-surface T,r=o and they are connected to the 
Lagrangian comoving coordinates a*(0, c?) at r = by the change of coordinates B i-^ 5(0, a). 
By inverting B = E(r, cTq) to Bg = g^{T,B), from Eq.(3.2) wc get a*(r, c?) = «*(0, ^g(r, a)). 
While on S,- with r > the position of the flux lines is identifled by the Eulerian coordinates 
B — E(t, Bo) in the Eulerian point of view, in the Lagrangian one this position is identifled 
by the Lagrangian coordinates q;*(t, B) — q;'(0, (t, B)) [see also Eq.(3.9)]. 

If we remember the deflnition (2.12), we have 

da%T,E{T,Bo)) da'{T,E{T,Bo)) jr/^<^^^:i\\M ^^^i^^^o) n Ci-i^ 
= q:^ + ^(^' '^o))\ f q:^ — = 0, (3.3) 

so that 
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= (r,E(r,..))]. = (3.4) 

Since we have 

inr, S(r, <7„))] , ^ . - (3.5) 



o o 



then we get 



J^(t, E(t, = - det (/(r, E(t, = -det"^ f |^ ) dct ' 



= Uoiao) det-' l^j . (3.6) 



The function 



n„(a,) = -det(^^^), (3.7) 

is the particle numerical density on the Cauchy hyper-surface I^t=o O'l^d is known from the 
fluid initial conditions. With the previous results, we can write 

.1 / <9S \ 



J'"(t, E(t, ao)) = no{ao) det" 



\d(To ' 



J^(r,S(r,a,)) =n„(a,)det-i (^) ^^'^^ 



Let us notice that, as said in the Introduction, any functional of the thermodynami- 
cal functions admits many expressions Q{T,a) = Q{z{T,a)) = ^(r, S(r, cTq)) = ^(r, cTo) = 
Qa{T,S{T,(T)) (we have added the expression in terms of the Lagrangian comoving coor- 
dinates) ^. In particular G{T,ao) and Q{t,B) are the expressions in the Lagrangian (or 
material) and Eulerian (or local) point of view respectively. While Q{t,Bo) is defined 
only on Et-=o, G{t,B) gives the expression on arbitrary hyper-surfaces S^^o- We have 
^('r>^)l?=s(r,ao) ^ ^{T,ao) and ^(r, ct) = g{T,ao)\ao=gs{T,a)- Therefore we get the follow- 
ing equation connecting them 



^Note that in general each of these functions is actually a functional of the associated coordinates 
and of their time and spatial gradients. 
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g{T, B) = j d^ao det S^{a - E(r, a<,))^(r, (3.9) 

The particle density fi is still given by Eq.(2.18) with the J"^ given by Eqs. (3.8). Then, 
introducing the notation 



T^i^, ^o) = y ^rr(T, S) + 2 C/^^(t, + frs(T, S) ls=E(r,a„)' (^-lO) 

we get 

n(T, ao) = n(r, a)|-^g(^^-^) = no(ao) ^ ^ = ^(^, ^o)- (3-11) 

^(r, E(r,ao)) 



Instead the entropy per particle s can be rewritten as a function dependent on &o alone 

s(q;(t, E(t, a^))) = s{ao{ao)) = So(ao), (3.12) 

and it is given with the initial conditions. In this case the constraint (2.9) is trivially satisfied 
being equivalent to 

|-s,(a„) = 0. (3.13) 
Analogously the particle number constraint is satisfied by construction being in particular 
M^l d^ar{r,a)^ [ d^aoUoiBo). (3.14) 

The action in the new coordinates is 
S^-Jdr d'a ^gir, a) 

^d'aoS'{a-E{T,ao)) det p{n{T, a^), s,{a,)) ^ 

= drd^ao\Jg{T,T.{T,ao)) det (^ "^^^^^"^ j p{Ht^^o),So{Bo)), (3.15) 

with the spatial integral restricted to the volume V^(0) in E^-^q- 

The equations of motion of the fluid variables derive from the stationarity of the new 
action with respect to the variations 5S(r, cTq). If the four- velocity f/^(r, ^o) is defined by 
Eqs. (2.10) with the new expression for the J^(t, Bo) given by Eqs. (3.8), that is 
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t/^(T,a„) = — i^^^^i^, where S^ = (r,SO, (3.16) 
then the equations of motion are 

^^'^"^ da'' Q:^^fi^^^o) + U {r,(To) -Q^V^{T,ao) -T{T,ao) ^0, (3.17) 

with V^{t, do) = flij, do) U^{t, (To). Let us prove that these equations are equivalent to Eqs. 
(2.19). First we can observe that from Eq.(3.9) we get 



da 



.g{T, d) = J d'ao det (^^4^ j ^'^^ " a^^^^' ^"^^ 



^g{r, a)^j d'ao det (^^^j 8\a - E(r, Bo)) x 

x(^-^^i^V,^^.)^l^(r,.^.), (3.18) 



o , 



where we have adopted the notation 

K:^r,ao)^-^^^^5t. (3.19) 

Using the rules (3.18) on the equations of motion (2.19) we can verify that Eqs. (3.17) 
are the material (Lagrangian) representation of the local (Eulerian) equations (2.19). ^ 
Therefore the action (3.15) defines the correct fluid equations of motion in the Lagrangian 
(material) point of view. 

For the stress-energy tensor we get 



^Let us observe that we get 



D ^, dg{T,a) d 



dfcFo det — — j 5 (o- - E(r, cTq)) 



dT ' 

with^;(T,a) = (^(r, 

i?o)) 1^0=92 with denoting the material temporal derivative. 
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I(r, 



= / d'oo det (^^^) - E(r, Bo)) f^^ir, a^), (3.20) 



where 

f^^{T,ao) = {p{T,ao) +p{r,ao)) U^{T,ao)U^{T,ao) - p{T,ao) g^^{T,^{T,ao)). (3.21) 

Obviously Eqs.(3.17), being equivalent to Eqs.(2.19), imply the conservation of the stress- 
energy tensor (3.20). Again there are no equations of motion for the z^{t, a)'s, so that they 
remain gauge variables also for the action (3.15). 
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IV. THE HAMILTONIAN FORMULATION. 



Let us study the Hamiltonian formulation of relativistic perfect fluid dynamics implied 
by the action (3.15). In this Section we use the notation /(S) instead of /(r, E(r, CTo)) for 
the sake of simplicity. The action (3.15) depends on the generahzed Eulerian coordinates 
E(t, a) of the fluid and on the embeddings z'^{t, a) as configurational variables. 

A. The First Class Constraints. 

We can define the following canonical momenta: the momentum density for the fluid 



Kf{T,ao) = 



5S 



n 



(4.1) 



and the momentum density of the embedding 

5S 



P/x(^, ^) 



+ 



52;^ (r, a) 

[ d^Uo det , ^ 
J \ auo 

j d^aoS{a -T,{T,ao))no{ao) 



dh J dz^{T.) 



+ 



dp z,,{T.) + z,,{T.)^ ^ 



dh n 

The Darboux canonical basis of phase space is represented in the following table 



(4.2) 





E^(r, 











(4.3) 



and we assume the following non null Poisson Brackets 



{17{t, ao),K,{T, a;)} = -51 5^3, - a'^). 



(4.4) 



As we said at the end of Sections II and III, we expect that the momentum density (4.2) 
is equivalent to four first class constraints. Since we have 
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z^.iY) = 0, 



then, using Eq.(4.1), Eq.(4.2) implies 



P(E) = V7(S), 



(4.5) 



p^(T,a)P(T,a) = j d'ao8\a-t{T,Bo)) x 



X 



det 



'9E(r^ 



7(S) + 



+ 



1 



^det(g)^9^ 



(^,0^0)- 



(4.6) 



The first expression is aheady in a constraint form. In the right-hand side of the second 
one the only dependence on the velocities is in the density n(T, Bo) [see the definition (3.11)]. 
Nevertheless we observe that from the definition (4.1) we obtain 

-f''^{T)Kf{T,ao)K's{T,ao) = \ no{Bo) ||(^'*)J 



nl{Bo) ^ 



+ 1 
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Then we can replace n(T, Bo) with the (implicit) solution X{t, Bo) of the equation 



(4.7) 



l'\T)Kr{T, Bo)K,(r, Bo) = (no(Bo) §j^{X(r, Bo), s) 



nliBo) 



X\r,Bo) l^det(g) ^)^ 



+ 1 



(4.8) 



It is evident by inspection that X{t, Bo) is a function of the canonical variables alone. In 
other words X{t, Bo) is independent from the r-derivative of the z'^'s and E's. It depends 
on the initial particle density 77.0(^0), on the Eulerian coordinates through — , and 
on the Eulerian momenta through 7'''*(E)K'r(r, Bo)Ks{t, Bo). 
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Then we have the following Dirac constraints: 
Hf{T, a) = Pi^{t, a)z^{T, a) - £ao ■ S^{a - S(r, ao))Kf{T, do) ^ 0, 



n^r, o) = p^(t, a)F{T, a) -J d'a, ■ 5\a - E(t, a,)) x 



det 



1 dp 



7 



(s)det(£)A-ax 



(r, ao)^0. 
(4.9) 



As shown in Appendix B, following Ref. [6], only in few cases, including the dust, the 
photon gas and a polytropic with n=|or7=l + ^ = 3, a closed form of Eq.(4.9) can be 
obtained. 

Working on the imphcit definition of X{t, Uo) given by Eq.(4.8) and using the results of 
Appendix C, it is possible to prove that the previous constraints satisfy the Dirac algebra 



{n^{r, a),n^{T, a')} = H,{t, a)^'\T, a) + Hr{T, a')Y\T, a') — 5\a - a'). 



d 



da' 



d 



{n^{T, a),nr{r, a')} = n±{r, ^')g^S%a - a'), 

so that they arc first class constraint. 

It is convenient to rewrite the constraints in the form 

W{t, a) = n±{T, a)P(r, a) + 7i,(r, a)Y'{T, a)z^{T, a) = 
= P^{t, a) - P(t, a) J d'ao ■ S'{a - E(t, a^)) x 

- z'^ir, a) I d^uo ■ 5\B - E(r, Bo))Kf{T, Bo) ~ 0, 
because then we get 



(4.10) 



(4.11) 
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{n^{T,a),n'{T,a')}^Q. 



(4.12) 



The Hamiltonian gauge transformations generated by these constraints change the form 
and the coordinatization of the space-hke hyper-surfaces E,-. Therefore the embeddings 
z^{t, a) are the gauge variables of this special relativistic general covariance according to 
which the description of isolated systems (here the perfect fluid) does not depend from the 
choice of the 3+1 splitting of Minkowski space-time. 

The Dirac Hamiltonian is Hd = J rf'^a A^(t, a) 7i'^(r, cr), where the A^'s are arbitrary 
Dirac multipliers. Since only the embedding carries Lorentz indices, the generators of the 
Poincare' group are Pfx — J d^a Pix{t, o) and Jf^^ — J d^a (2;'^ p'^ — p^^) (r, a). 



B. The Restriction to Space-Like Hyper-Planes. 

Following Refs. [8,10,3] let us restrict ourselves to foliations whose leaves are space-hke 
hyper-planes by adding the gauge-fixings: 

C^{t, a) = z^{t, a) - x^ir) - 6^(r)(7' ^ (4.13) 

With this condition many geometrical quantities take a trivial expression, in particular 
we have: 

z^{T,a)^b'^{T), 
z^{T,a)^x''{T)+b'^{T)a^, 

9fs{r,a) ^ -Sfs, Y'{r,a) ^ -6''\ 7(r,a)Ril. (4.14) 

We also introduce the following notation for the unit normal to the hyper-planes defined 
in Eq.(A3) 

6^!(t) = P(r) P(t, a) (4.15) 

The hyper-planes define a true global foliation only if the normal is r-independent, 
because only in this case the hyper-planes are parallel and not intersecting and there is a one 
to one global correspondences between points z^ and coordinates r, a. We ignore from now 
on this observation, but we shall return on the consequences of the P = constant request 
at the end of this Section. 
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Since we have 



{C'^(r,a),H.(r,a')} = -rj^J\a - a'), (4.16) 
it is possible to define the foUowing Dirac brackets 

{A, BY = {A, B}- j d'a ■ [{A, C{t, a)}{n^{T, a),B}- {A, H^ir, a)}{C,iT, a), B}] . 

(4.17) 

In the reduced phase space the hyper-surface canonical variables z^{t^(j)^ p^{t,(t) are 
reduced to only ten canonical pairs and the first class constraints (4.11) are reduced to only 
ten [8,6]. 

We can verify that four pairs of canonical variables are the centroid x'^{t), used as origin 
of the 3-coordinates on the hyper-planes, and the conjugate momentum p'^(t) 

{x^{T),p''iT)r^-rj^r (4.18) 

Since p'^ is the Poincare' generator of the translations, it describes the total 4-momentum 
of the system. As a consequence we have the following decomposition for the canonical 
generators of Lorentz transformations inside the Poincare' algebra 

j/^'^(r) = x''{t)p%t) - x%t)p''{t) + S^'^ir). (4.19) 

The remaining canonical variables defining the hyper-planes are the variables 0a (t) ( 
A = 1, 6) that parametrize the orthonormal tetrad b'^{T) such that 

&!^(T)r7;..6^(r) = vAb: (4-20) 

and the associated conjugate variables 7a (r). Nevertheless it is possible [8] to use a set of 
redundant, non independent and non canonical variables. These are the tetrads &^(t) and 
the spin tensor S'^^{t) of Eq.(4.19), if they satisfy the following Poisson brackets ^ 

{^'^'^(t), SP^iT)}* = CZ'^S'^^{t), (4.21) 



They are the Dirac brackets enforcing the fulfillment [24] of Eqs.(4.20). 
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with the Lorentz group constant structures: 

f^ixvpa _ v p na _ p p va _ u a pp , p a up 

Finally we have the unchanged fluids variables: 



(4.22) 



(4.23) 



In conclusion this non-Darboux canonical basis of the reduced phase space can be rep- 
resented in the table: 





E^(t, 


Bo) 






Bo) 



(4.24) 



The stability in time of the gauge fixings (4.13) force the arbitrary Dirac multiplier to 
take the reduced form: 

d 



dr 



((r, B)^0^ X^{r, B) ^ A''(t) + X^''{r)r),pb',{r)a' . 



(4.25) 



The Dirac Hamiltonian become 

Hn = X^{r)H^{r) + \^''{t)H,,{t), 

where, if we define 



(4.26) 



A.(.)^/.3..|det(g 
V'{t)^ I d^aoK'{T,Bo), 
j-(r) = ]C'{t) = - I d^ao^'{T,Bo) x 



{r,Bo), 



X 



det (g) X ax 



{r,Bo), 



J'Ht) 



) = j d?ao [Tr{T,Bo)K\T,Bo)-T.\T,Bo)K'{T,Bo) 



(4.27) 



we have the following form of the ten remaining first class constraints 
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//^(r) = p/'ir) - F{t)M{t) + }ff{r)V'{T) ^ 0, 

= S^'{t) + [6^(T)r(r) - 6,^(r)r(r)] J^'\t) + 

- ^ mrmr) - 6^(r)6^(r)] J^^(r) 0. (4.28) 
They satisfy the Poisson algebra 

{H>^{t),H'^{t)} = {H>^{t),H-^{t)} = 0, 

{//^'^(r), Hf'ir)} = C^J^mf^ir) ^ 0. (4.29) 

These ten first class constraints imply that x^{t), bf{T) are gauge variables and that the 
description is independent from the choice of the space-like hyper-planes. 

We can discuss now the condition = constant. 

By using the brackets (4.21), the Hamiltonian (4.20) implies the following equations of 
motion for the normal /^(t) 

^Ur) = 2A^4r)/'^(r). (4.30) 

Then the condition If^ — constant restricts the arbitrariness of the A^,/(t)'s with the 
condition 

A^.(T)r = 0, ^A'^'^(t) = 6^(t)6,^(t)A-(t), A-(t) = -A-(t). (4.31) 

In other words the condition = constant have to be interpreted as 3 gauge fixing condi- 
tions for the Lorentz transformations generated by the constraints H^'^{t) 0, leaving only 
the possibihty of Hamiltonian gauge rotations generated by Hfs{T) — b^{T) b'^{T) Hn,y{T). 
This implies that a fohation with space-like hyper-planes has 7 gauge variables: x^{t) and 3 
angles inside 6f! (r) describing the linear acceleration of the origin and the arbitrary rotation 
of the spatial axes with respect to an inertial frame, respectively. Correspondingly only 7 of 
the first class constraints (4.28) are independent, namely H^{r) and Hfsir). 



26 



V. THE REST FRAME INSTANT FORM. 

In this Section, following Ref. [8], we define a new instant form of dynamics [24] called 
the rest frame instant form. This is done by selecting all the configurations of the isolated 
system with time-like conserved total 4-momentum and, for each of them, by choosing 
the foliation whose space-like hyper-planes are orthogonal to this 4-momentum. Physically 
these hyper-planes, named Wigner hyper-planes, correspond to the intrinsic rest frame of the 
configuration of the isolated system. Moreover we have to fix the four acceleration degrees 
of freedom of the centroid a;^(r), reducing it to the world-line of an inertial observer. 

Before doing this we have to recall the notion of standard Wigner boost and Wigner 
rotation. It is known that in the rest frame a timelike four- vector assume the standard 
form p'* = (Vp^, 0). In the general theory of the induced representation of the Poincare 
group a standard Wigner boost L^^^{p,p) for time- like Poincare' orbits is defined such that 

L>',ip,p)p''^p'. (5.1) 

We can use the rows of the standard Wigner boost matrix to define a orthonormal tetrad 
(we define n^^ — p^/ Vp^) 

e^(p) = L'^„(p,p) = n^ 

e^,{p) ^ L\{p,p) = (-n,;5;-nV(l+no)^/'). (5.2) 
We use the notation 

The inverse of standard boost defines naturally the inverse tetrad 
such that 




elip) = {d^'us, 51 - 5''5jnn'n,{l + n„)V^), (5.5) 

and 
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(5.3) 



(5.4) 



r 

Moreover we can verify that 



(5.6) 



p'i-Ap)=p'i-,^AiP)=o. 



(5.7) 



The standard boost can be used to define the Wigner rotation R{A,p) associated to any 
Lorentz transformation A 



^ ^ 



L \p,p)A ^L{Ap,p) 



R{A,p) 



(5.8) 



Then we have the following property 

e(f(Ap)=AVr(p)i?%(A,p)- 



(5.9) 



We can use the tetrad (5.3) to define a gauge fixing for the &^(t). Before doing this, it 
is useful to replace the centroid x^{t) with the new one 



and the table (4.24) with the following set of non canonical variables 







S^(T,ao) 






P^{t) y\{T) 






The q^{T),p''{T) and E^(t, 




are still canonical 




(t)} = 0, 




r)} = 


{?"(r),S^(r, 


^o)} = 0, 





(5.10) 



(5.11) 



(5.12) 



but for S^^{t), &i(T) we have the following non canonical brackets 



28 



{g^(T), S'^^ir)} ^ 0, {q'^ir), ^^(t)} ^ 0. (5.13) 

Eq.(5.10) also implies the following new decomposition for the canonical generators of 
Lorentz transformations 

JM'^(r) = g^(T)p'^(r) - q''{T)p^{T) + VI^''{t), (5.14) 

where 

Si-(r) = 5- - i.„-(p),.B (^p- - ^p-) S"". (5.15) 

Since we have 

{g^, i^"} = 0, {g^, n"''} ^ 0, (5.16) 
we see that is not a true four- vector. 

If we restrict ourselves to configurations with p^{t) time-like (p^(r) > 0), there exists a 
family of space- like hyper-planes orthogonal to p'*(t) and we can select this family with the 
gauge fixing [8] 

n{r) = 6La(^) - ^a(p) ^ 0, (5.17) 

where the index f is enforced to coincide with r with transformation property given by 
Eq.(5.9). These gauge fixings imply A^,y(T) f« ^ and we have the Dirac Hamiltonian 

Hd = \^{r)H^{T), (5.18) 

where 

H'^ij) ^jf- n'^M(T) + e(f(p)7'"(T) ^ 0, (5.19) 

are the remaining four first class Dirac constraints saying that x^{t) and terefore g^(r) are 
gauge variables. Since we have x^{t) = {x^{t), Hd} — —X^{t), we see that the centroid has 
an arbitrary (gauge) acceleration described by the remaining Dirac multipliers. It is useful 
to observe that we can rewrite the constraints and the Dirac Hamiltonian in the form (we 
choose the positive sheet of the mass hyperboloid) 



^So that Eqs.(4.31) are satisfied. 
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p-^lr) - 6;;(p)i^'^(T) « 0, 



= A(r) 7i - A(r) • V, 



where, recalling the definitions (4.27), the quantity 



(5.20) 



n^(CTo) 1 dp 



X dX 



is the invariant mass and where 

^(t) = I dV<,X(T,a<,)«iO, 



(5.21) 



(5.22) 



is a constraint on the total momentum of the fluid inside the Wigner hyperplane: it gives 
the rest-frame condition. 

See Appendix B for the expression of the invariant mass M. for the dust, the photon gas 
and the polytropic with n — \. 



After this gauge fixing the non canonical variables 5''*'' can be everywhere substituted 
with the expression implied for them by the constraints H^^{t) a; of Eqs.(4.28). By 
defining 



^-(r) = j d'ao [E^(r, ao)K^{r, Bo) - E^(r, BoW{t, Bo)] , 
j-(r) ^JC'^-J d^ao S'-(t, Bo) x 



det 



dan 



7(S) ip{X,s)- 



{r, Bo) 



^-J d'ao^'{T,Bo) A{t,Bo), 



(5.23) 



we get 
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(5.24) 



After the elimination of the variables 6^, 5''*'^ with Eqs. (4.28) and (5.17), the reduced 
phase space is spanned by the variables 





E'-(t, 











(5.25) 



Its canonical structure is defined by the new Dirac brackets 
{A, BY* = {A, BY - [{A, H^^YKp(p){Ta, BY + {A, nYKMiH"', BY] , (5.26) 
where 

KM = liv^p^ip) - Vu^ip)]. (5.27) 
The variables (5.25) are canonical with respect to this bracket 

{E'-(t, Bo), K^{r, a'JY* = <5(a, - a',), 
{q^, q'^Y* = 0, 

{q^,E^{T,ao)Y* = 0, 

{q^,K%T,ao)Y*^0. (5.28) 

On the contrary, with the new brackets the old centroid is not canonical, because we 
have 



{x^ E'-(t, ao)Y* + 0, {x\ K\t, a„)}** + 0, 

and this explains why we introduced the centroid q^ij). 
After the gauge fixing we have 



(5.29) 



and the canonical generators of Lorentz transformations become (v^ ~ A1) 



(5.30) 
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J''\t). (5.31) 

This form of the canonical generator of Lorentz transformations tell us that on a function 
dependent on the fluid variables alone F(E(t, ^o), i^(r, ao)) a Lorentz transformation acts 
as a rotation inside the hyperplane. This rotations is the Wigner rotation associated to the 
infinitesimal Lorentz transformation 

5F = d(jj^,{F, J""}** = dip''{F, J-}**, (5.32) 

where 

^o;^!/. (5.33) 



+ 



p° + Vp^ 



The canonical generators of Poincare group p'^{t) and J^''{t) are called external and 
they realize the true Poincare symmetry in the reduced phase space. Because of property 
(5.32), the Lorentz covariance of the theory is replaced with the Wigner covariance of the 

3- dimensional instant form variables on the hyper-planes. In this sense the only canonical 
variable with a non covariant transformation property is the pseudo four- vector q^. As 
said in Ref. [14], the centroid q^^ is interpreted as a non- covariant, hut canonical external 

4- center of mass. In particular its spatial components arc proportional to a 3-ccntcr-of- 
mass-likc position that is the classical analogue of the Newton- Wigner position operator, 
whose reduced covariance corresponds to the little group of time-like Poincare' orbits. 

For these particular properties the hyper-planes defined by the condition (5.17) have 
been called Wigner hyper-planes in [8] . 

We can also construct another realization of the Poincare Lie algebra. This reahzation, 
called internal, is constructed using only the fluid variables S(t, CTq), X(r, (Tq) living inside 
the Wigner hyperplane. In fact, from Eqs. (5.21), (5.22) and (5.23) the 10 functions: 
P,M,J,IC, where 

jr ^ ^^ruvjuv^ ^ jrr ^ (5 34) 

are the generators of the following canonical realization of the Poincare Lie algebra 
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{M,V'} = [VW] = 0. 



(5.35) 



This internal realization is unfaithful due to the constraints P ~ 0. These constraints 
says that three degrees of freedom, playing the role of an internal 3-center of mass of the 
fluid inside the Wigncr hyperplane, are gauge variables. As shown in Ref. [14], /C ?a are 
the natural gauge fixings to eliminate the internal 3-centcr of mass and to imply X[r) — in 
the Dirac Hamiltonian of Eq.(5.20). Only the invariant mass Ai and the rotation canonical 
generator are non vanishing and they appear in the external Poincare generators (5.31). 

With the gauge fixings /C 0, the Dirac Hamiltonian is reduced to Hd — A(t) H. 

To complete the definition of the rest-frame instant form we identify the temporal pa- 
rameter T with the common Lorentz scalar rest frame time of the centroids q'^ and x^, 

T{t) = n^x^ir) = n^q^ir). This is realized by introducing the gauge fixing 



T{t) ^ T, 

implying A(r) = —1. By using the canonical transformation 



(5.36) 



z{t) = Jp^ ■ Q{t) = Jp^ 



with inverse 

g°(T) = + P I T{t) + 



k ■ z(t) \ . z (, , k ■ z(t) , , 



(5.37) 



p" = S\/l + k^, p = £k, 



(5.38) 
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we arrive at the following Darboux canonical basis 
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T{r) 


z{r) 


S{t) 





(5.39) 



In the rest frame instant form we have T ^ r, £ ~ Ai, Ho = and the external 

canonical non-covariant 4-center of mass q'^ is interpreted as a decoupled point particle clock. 
However, since the gauge fixing T — r ^ is explicitly r-dcpcndcnt wc get that the effective 
Hamiltonian to reproduce the Hamilton equations for E(r, ao), K{t, do) when A(r) = —1 is 

H = M- A(r) • V. (5.40) 

In the new canonical variables (5.37) we have [Q* = {l/2)e^^^Q,^^ with Q,^^ given by (5.30)] 

J(r) = z{t) X k{T) + Q(t) = Q(t) X p{t) + Q(t). (5.41) 

The reduced phase space now is 



z{t) 


E(r, do) 


k{T) 


K{r,do) 



(5.42) 



The final result is a new instant form of dynamic [25] , the Wigner-covariant rest frame 
instant form. In this form the system is described by a reduced phase space (5.42) formed 
by two sectors: 

i) the sector of the external decoupled point particle clock described by the canonical 
non-covariant variables ^(r),A;(r); 

ii) the sector of the internal Wigner-covariant variables T,{T,do),K{T,do) (they are 
Wigner spin-1 3- vectors), living inside the Wigner hyper- planes. Since they are restricted 

— * — * 

by the constraints V ^ 0, K, 0, only variables relative to an internal inessential 3-center 
of mass (see next Section) are physical. 



= z/£ is the classical analogue of the non-covariant Newton- Wigner position operator. How- 
ever we cannot replace 2 and k with Q, p, because they do not have vanishing Poisson bracket with 

the internal variables S(r, cJo) and K{T,do). 
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The external variables define the 4- unit vector — {\ 1 + k^,k), which identifies the 
direction of the total four-momentum with respect to an external inertial observer. 

In the rest-frame instant form of dynamic the 10 canonical generator of the external 
Poincare group are 



k xQ 



K ^ -zyi + k'^ - 



1 + Vl + P 



^ -Q - -—^^^==, (5.43) 



where 



.2 

is the (interaction-free) spin with respect to the external center of mass. 

The properties of the instant form are now shown explicitly: i) p, J do not depend on 
the interactions; ii) only the 4 generators p°, K depend on the dynamics through Ai. Even 
if in a generic instant form the dynamics is determined by four independent potentials (the 
dynamical SU (2) of Ref. [26]), in the rest frame instant form there is a unique function, the 
invariant mass Al, carrying the whole dynamical information. 

Let us observe that the ea;terna/ boost generator K'^ = of Eqs.(5.31) can be rewritten 
either in the form 

K{t) = g^(r)p - q{T)p'^ - (5.45) 

or in the form 

K{r) = q\T)p + K'{t). (5.46) 
If Q{r) is defined by (5.37), we have 

^^^^ = - Ti^wwy 
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and 

Hr)=-^-4^^. (5.48) 

In this way we recover the usual definitions of the canonical, non-covariant relativistic 

3- center of mass [27]. The first is given using the complete external boost generator, the 
second using the boost generator on the hyperplane ^"(t) = 0. Being n^x^{T) = n^j^q^ir), the 

4- center of mass ^''(t) defines a point on the Wigner hyperplane different from the centroid 
x^{t) — z'^{r,6), with coordinates a{q) ^ such that q^ij) — z^{r,a{q)). Using only the 
canonical generators of the external realization of the Poincare group it is also possible to 
define the external non-covariant, non-canonical M0ller 3-center of energy [27]: 

Sir) = (5.49) 

and the external covariant, non-canonical relativistic Fokker-Pryce 3-center of inertia [18]: 

#M_Pi^. (5.50) 

In Ref. [14] it is shown how to identify the position on the Wigner hyperplane of the 
associated external (pseudo-vector) 4-center of energy and of the external 4-center of 
inertia V^, which is a 4- vector by construction. If we put q^ = {q°]q), then we have 
j^n _ i^qO-^R-\- qo and — {q°]Y -\-q° p). In Ref. [14] it is also shown that all the possible 
pseudo- vectors qi^ and fill a world-tube around the 4- vector y*, whose radius is the 
M0ller radius [28] p = |^ of the fiuid configuration. This radius is defined by the Poincare' 
Casimirs of the external Poincare' group (p^ = Al^, — —M^O,'^) and is a classical unit 
of length determined by the Cauchy data of the configuration of the system. See Ref. [3] for 
the properties of this radius and for the proposal of using it as a natural physical ultraviolet 
cutoff at the quantum level for all rotating configurations of the isolated system. 

Let us conclude this Section with the Hamilton equations associated with the Hamilto- 

— * 

nian (5.40) in a gauge where A(t) = 

= {S(r, cro),M}, 

^^^^^{K{T,ao),M}. (5.51) 
For the dust Eq.(B6) implies 
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E(r, ao) = a. + , ^"l r. (5.52) 

Let us remark that in this description these hyperbohc Hamilton equations replace the 
hydrodynamical Euler equations (2.26) imphed by the conservation (2.23) of the stress- 
energy tensor. 



37 



VI. INTERNAL CENTERS OF MASS AND RELATIVE VARIABLES 



In the rest frame instant form defined in the previous Section the 3-dimensional variables 



To select the physical degree of freedom is a problem equivalent to determine the internal 
3-center of mass and relative variables on the Wigner hyperplane. As already said, the 
internal 3-center of mass variable on the Wigner hyperplane is a gauge variable, because 
the role of true center of mass is played by the external, canonical non-covariant, 3-center of 
mass z. On the contrary the relative variables will be the physical variables and they will 
describe the reduced phase space after the final gauge fixings, whose natural form is /C ~ 0. 

To justify the gauge fixing /C ~ we must perform two steps. First we select a naive 
internal center of mass position X canonical with respect to the total momentum V and the 

— * 

associated canonical relative variables. The naive center of mass position X allows to define 

— * 

the gauge fixing A" fa 0, but this gauge fixing has the unpleasant property that the arbitrary 

— * — * 

Dirac multiplier A(t) is fixed to a non null value (A(t) ^ 0). Then we use the internal 
relative variables obtained as auxiliary variables in the first step for defining the relative 
variable respect to a 3-center of mass Q such that the gauge fixings Q ~ (identification 
of the internal 3-center of mass with the centroid x^{t) origin of the 3-coordinates) imply 
A(r) ^ 0. This is done using the method of the Garien/iaits-^'c/iuiarfo transformation of Ref. 
[29] (see also Appendix D). 

Let the total internal 3-momentum and an internal naive 3-center-of-mass position on 

— * 

the Wigner hyperplane X be defined as 



— * 

on the Wigner hyper-planes are not all physical due to the first class Dirac constraint ~ 0. 





(6.1) 



with 




(6.2) 



and 



(6.3) 
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Let the internal relative canonical variables 3?^(t, CTq), P^{t, Bq) be defined in such a way 
that we get 

K'{T,a,) = ^^V^{t) + J rfV>'-(r,a;)r^K,a,). (6.4) 

The kernels Tk, will be specified by some conditions that we will analyze shortly. 
From the definitions (6.1) we obtain that 

j d^aono{ao)TY:{ao,a'^) = 0, 

J d^aorK{a'„ao)^0. (6.5) 

We impose the following canonical property 

{W{t, ao), p^T, = S^'S{ao - a;). (6.6) 
By using Eq.(6.3), we can verify that from Eq.(6.4) we obtain the canonical property 

{E'-(t, Bo), K^ir, <)} = 5- 5{ao - (6.7) 

if 

J (fao Tj:{aoi, a)rK{a, ao2) = - ^"^"^^ + S{aoi - ao2)- (6.8) 

When Eqs.(6.8) hold, Eq.(6.4) is consistent with the following definitions 

r(r,ao) = J rfV;rx(ao,a;)S'-(r,a;), 

p^{r,ao) = / ciV;X^(T,OrE(5^;,?o). (6.9) 

In fact, if we substitute these expression in Eq.(6.4), using Eq.(6.8) we have an identity. 
Again from Eq.(6.9) and using Eq.(6.7), we can find Eq.(6.6) if the following condition is 
verified 

j d^ao TKiBoi, ao)Ts{ao, ^^2) = ^(^01 - ao2)- (6.10) 
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Eqs. (6. 5), (6. 8) and (6.10) are a set of conditions that have to be satisfied by the kernels 
r. These conditions are not independent: it can be proved that Eqs. (6.8) e (6.10) imply 
Eq.(6.5). 

We can also verify that, using Eq.(6.5), we have 

J{t) = X{t) X V{t) + 1 Suo 5ft(T, Go) X p(t, Go), 

)C{t) = -M{t)X{t)- J d'God^G',rj:{Go,Gj^{T,G'jA{T,Go). (6.11) 

Let $n(ao) be a base of orthonormal functions on with n = (ni,n2,n3) a set of 
multindices. Then we can consider the coefficients 

rn{r)^ J d'Go^n{^o)^{r,ao), 

Pnir) ^ J d^Go^n{ao)p{r,Go), (6.12) 

such that 

^(r,Go)^J2^n(r)^n{ao), 

n 

P{r, ^o) = Y.Pn!yr) *il(^o). (6.13) 

n 

Moreover from Eq.(6.11) we get 

I d^Go ^(r, Go) X p(t, Go) = ^ 4(r) X p„(r), (6.14) 

and 

K(r),p^(r)} = 5-5^. (6.15) 

In conclusion the coefficients fn{T), Pnij) are a set of infinite canonical variables that 
we can use as internal relative variables. These variables are useful for defining the canon- 
ical transformation that will realize the separation between rotational and shape degree of 
freedom in the next Section. For the time being we can use them to rewrite the definitions 
(6.4), (6.9) in the form 
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E(r,a,) = ;?(r)+^rJ(a„)4(r), 

n 

rn(T) = J dV„rJ(a„)E(r,a„), 

Pn(T) = y (iV„rJ(a<,)^(T,a„), (6.16) 



where 



rj(a„) = / dV;$^(a;)r^(a;,a,). (6.17) 

The conditions (6.5) and (6.8), (6.10) on the kernel functions F are rewritten in the form 
I £ao rioiao) rj(a„) =0, j (a„) = 0, 

X] rJ(CTio)r^(CT2o) = + (5(5^^^ _ 5^2^)^ 

y d?ao rj(ao) r^(^o) = <^nm- (6.18) 

Some possible solutions for the kernels V are derived in Appendix E. 

As said at the beginning of this Section, the internal 3-center of mass-like position X is 
such that the gauge fixing X Q does not imply A(t) as can be checked by using the 

— * — * 

Dirac Hamiltonian Hd — Ai — A(t) • V in the gauge T t. We want to replace it with 

— ♦ — ♦ -* 

another internal 3-center of mass Q such that the conditions imply A(t) ^ 0. 

To this end we construct the internal 3-centers of mass, energy and inertia in analogy 
to the external ones of the previous Section, using the internal realization of the Poincare 
algebra instead of the external one 
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CM = .m — ^"xf» . « 

W = - ^i±L . (6.19) 



where 



J(r) = Qx:P + 5(r), 

/C(r) = -M(t) ;e(T) - / rj(a„) A(t, a,) = 

= -m{t) n{T) ^ -m{t) q(t), 

'^(T)=^(r)-i:^/^'^orJ(c?„)A(r,cx.). (6.20) 

These three centers are weakly equal due to the constraint P pa and they are all 
canonically conjugate to V 

and such that 

{7^^ = {Q\ MY* = {y\ MY* = ^ ~ o. (6.22) 

— * 

But Q is the only one such that 

{Q'-, Q^Y* = 0, (6.23) 

namely it is the real internal canonical 3-center of mass. 
If we adopt the gauge fixings 

q(t) ^ n{T) ^ y{T) ^ 0, (6.24) 
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then we get 



Q(t) = {Q{t),M - A(r) • P} = -A(r) « 0, 



r„\T 



M(t) 



E(t, CTo) J2 



/rfV:r^(a:)A(r,a:) 



M{r) 



(6.25) 



Eq.(6.19) shows that Q(r) ^ is equivalent to the condition JC ^ 0. After this gauge 
fixing, in the internal unfaithful realization of the Poincare group there are only four non 
null functions: and Jj. 

In Appendix D it is show that with the relative variables 3?^(t, CTo), p*(r, a,,) we can 
realize the Garthenaus- Schwartz canonical transformation 







Q{r) 




E(t, ao) 


— > 








V{r) 





(6.26) 



In Appendix D it is also shown that, after having gone to Dirac brackets with respect to 

— * — * 

the second class constraints V ^ 0, Q fa 0, we get 



P=Q = ^'■{T,ao) = U'-{r,ao), p'^(t, a„) = p^(t, a„). 
Then the final reduced phase space is 



(6.27) 



z{t) 




kir) 





(6.28) 



where 



J{t) fti 5(r) = / ^(r, a^) x p(r, a^) = ^ r;(r) x p„(r). 



(6.29) 



On the Wigner hyper-planes the kinetic term K'^{t, do) appearing in the solution X(r, do) 
of Eq.(4.28) is 
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X^(r,a,) ^ Y: r^^(^o)rJ(a„)p^^(r) .p;^(r), (6.30) 
while the dependence of X{t, Bo) on the generahzed Eulerian coordinates is concentrated in 
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VII. ROTATIONAL KINEMATICS 



The internal relative canonical variables fn{T), Pnij) are Wigner spin 1 vectors under 
rotations. Then we can do on them a canonical transformation that generalize the results 
obtained in the N particles case in Refs. [13,14]. The r„ are interpreted as a set of infinite 
relative position vectors; we want use them to construct a dynamical body frame as done in 
Refs. [13,14]. For this we have to select a pair of these vectors. We assume to choose the 
vector positions with multindices UxiUa as preferred vectors 

f„. = Ri, pu. = Hi, i = 1, 2, (7.1) 

and we use them to define the orthogonal vectors 

-± Ri + R2 ^ Ri- R2 -t ^ , s 

N = , X = ^ , N-x = 0- (7.2) 

The most convenient choice of these two vectors will be dictated by the spatial form of 
the initial density no((?o)- 

Then a dynamical body frame is defined by the associated unit vectors and their orthog- 
onal complement 

W = (x(t),^(t)xx(t),7V(t)). (7.3) 

By construction this frame is a orthonormal frame that rotates with the motion of the 
fiuid. In this sense it generalize the concept of body frame of a non relativistic rigid body 
and we can apply the same observations and interpretation done in Refs. [13,14]. Moreover 
we observe that 

{7V^ N^} = {x^ X^} = {^^^ X^l = O. (7.4) 

All vectors can be projected on this dynamical body frame so to obtain their components 
on it. In particular, for the relative angular momentum, given by Eq. (6.29), its components 
on the dynamical body frame are 

r{T)=S{T)-br{T), (7.5) 

or more explicitly 

S'^S-x; S'^S-iNxx); S'^S-N. (7.6) 
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Using the results of Refs. [13,14] we can construct the following quantities 

vP = ^1 X rii - ^2 X lis, 



(7.7) 



Ri —\Ri\, Hj — Ilj • Ri, i — 1,2, 



(7.8) 



and, for n^Ui,U2 



rl = rn-x, rl = fn- N XX, fl = fn-N, 



Pn^Pn-X, Pl^Pn- N XX, pl^Pn' N. 

Then the transformation represented in the table 





> 


1 ^ 1 






1 N 1 


Ri 


R2 




Pn 




a 




7 




Hi 
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(7.9) 



(7.10) 



is a canonical transformation (the canonical pairs are the variables on the same column) if 
we define 

.1 {S X Nf 



a — tan 



(5 — tan 



[5 X (5 X 7V)]3' 
51' 



7 = tan =1", 
S 



W-{Nxx) 



(7.11) 



The canonical variables in the final basis have been separated in three sectors. The 
second and the third sector in the previous table are constituted by canonical variables 

scalar under rotations: these variables describe the shape of the fluid and we call them shape 
(or vibrational) variables (see Ref. [30] for the original definition of shape variables). The 
first sector is that of the rotational (or orientational) variables; these variables describe the 
rotational motion of the dynamical body frame. 
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It is useful to analyze the variables using a new orthonormal base: the spin basis of Refs. 
[31,13,14]. This basis is defined observing that there is only a unit vector ^ on the same 

— * — * — * 

plane of N and S, orthogonal to S such that 



a — tan" 



.1 {s X 
[Sx{Sx n)Y 



(7.12) 



Then the three unit vectors {IZ, S,S x TZ) are a orthonormal basis, the spin basis. By 
construction their components are given by the following relations 



' ^1 ^ y/S^-{SS)^ 



\s\ 



COS (3, 



C3 _ S-' 



(7.13) 



TZ = sin P sin a — cos /? cos a, 



7^2 = 



cos Psina — sin P cos a, 



^3 _ 



COS a. 



(7.14) 



{S xTZy = sin P sin a + cos /3 sin a. 



((S X 7^)^ = — cos P cos q; + ^ sin /3 sin a. 



(5x7^)3=^^^^sina. 



(7.15) 



We also define the angle ip such that 



cosip — — smip — — 



S 



S 



(7.16) 
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By definition of 1Z we have 

TV = 008-0 «S + sin 7t. (7.17) 

Moreover tlie definition of 7 implies tliat we liave 

S ■ x = sin i}: cos 7, 
5 • (x X TV) = sin'0sin7. (7-18) 

We complete the conditions on x, TV x x using the fact that the dynamical body frame 
and the spin basis are connected by a (proper) rotation 

X = sin tp cos 7 (S — cos ip cos 7 7?. + sin 7 5 x 7^, 
X X iV = sin'0sin7 «S — cos -0 sin 7 TZ — cos 7 S xTZ. (7-19) 

Substituting in Eqs.(7.17),(7.19) the expression given by Eqs.(7.13), (7.14), (7.15), (7.16) 
we obtain the elements of the dynamical body frame expressed as functions of the rotational 
variables alone. We can also define the Euler's angles of the dynamical body frame s^,e'^, 

cos£^ = N^; cose^ = , ^ ; cos£^ = 1 ^ ■ (7-20) 

They are as functions of the rotational variables alone. 

The Euler's angles (7.20) together with the relative angular momentum components 

5' = V^^^^^cos7, 

< 5' = V52-(5Vsin7, (7-21) 

define a non canonical transformation for the rotational sector. The corresponding canonical 
transformation is obtained using the canonical momenta 



As in Refs. [13,14] we adopt the y-convention of Ref. [32] 
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Pi — — sin £^ cos + sin £^ sin + cos £^ S^, 



P2 = sm e b + cos £ o , ps = o , 



(7.22) 



The inverses of the previous equations are 



^ cos s 

S =sin£^p2 r-pi + COS£^ COt£^P3, 

sin£^ 



2 sill £ 

S — COS p2 H 5- pi — sin £^ cot £^ p3, 

sin£^ 



-^3 



(7.23) 



and we get 



03 -23 "ol I • 2 • 3 "o2 2 "^3 

o = pi = — sm e COS £ o + sm e sm £ o + cos £ S , 



/ 1 7r\ sin£^ (— sin£^cos£^<S^ + sin£^ sin£^<S^ + cos£^<S^) — 

p = ( £ — — 1 — arctan =j =2 ^ 

2 / sin £2 (sin £^ 5 + cos S ) 



q; = arctan 



sin £^ (sin £^ 5^ + cos £^5^) \ S 



cos £^ I «S P — (— sin £^ cos £^ ^S"*^ + sin e"^ sin £3 + cos £2 <S^) 



=2: (5-/. 

This chain of transformations can be represented with the table 



(7.24) 



1 ^ 1 
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— > 






5^ 


)• 


Pi P2 P3 


a 




7 




£l 


£2 


£3 







(7.25) 



and with the following Poisson brackets for the non-canonical variables £*', [29] (/, g are 
functions only of these variables) 
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Of dg df dg 



^ / (9/ dg 



(7.26) 



where are the components of the right invariant vector fields on the group man- 

ifold of 5*0(3). The components of the dual right invariant one- forms are 



^ — sin £^ cos £^ sin £^ ^ 



2 ^3 cOS£^ 



V 



sm£ sm£ 

COS£^ 







-1 



if the Euler angles are defined by the convention R{e''') = R-i{s^) -R2(£^) R?,{i^^)- 

These transformations stress the canonical equivalence between the rotational variables 
and the canonical phase space of a non relativistic rigid body [32,33]; in other words, the 
rotational variables are the rigid body-like variables, whereas the shape variables describes 
the non-rigidity of the system. 

Using the results obtained in the three body case in Refs. [13,14] it is easy to construct 
the inverse canonical transformation and to express the original relative variables in terms 
of the rotational and shape variables. In particular it is trivial to observe that by using the 
rotation Rg{e^, s^, e^), for n 7^ Mi,2i2 we obtain immediately 



pl = R^,{e\e\e')pl 



(7.27) 



For 21 = Ml, 1^2 we can to use the three body results of Refs. [13,14]. In particular we 
have 



= 4, • (X X TV) = 0, 
rt. = 4^ • iV = i?, I iV I . 



(7.28) 



Then we obtain 



Tl^R^,{e\e\e')Tl, 



(7.29) 
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Finally, if we define 



we have 



I N I 
'2R~ 



5(12) + (-)^^'eVi 



Pi = (-) 



2Ri 



^ ^^+1 *^(12) ^ '5(12) 



p^^. = n, I TV I -(-) 



'5(12) + 



and then 



Therefore we get 
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VIII. THE INVARIANT MASS AND THE EQUATIONS OF MOTION. 



As we have seen, the invariant mass (Dixon's mass monopole as we shall see in the next 
Section) M. — J dPuo A(r, Gq) is the Hamiltonian on the Wigner hyper-planes in the gauge 

— * — * — * 

T T and Q ~ 0, where E(t, CTo) and K{t, Bo) are given by Eqs.(7.33). Therefore we have 
-^M — but it can be shown that ■^A{T,ao) ^ 0. Moreover, since /i„ in Eqs.(7.33) 
depends on both configuration and momentum shape variables, it can be shown that we 
have ^ 7^ 0, namely 

^5:^ = (<3(r)xi](r,a„))' + 

+ fl'.(.«(r)) E { ^ \vl(a,) - K[r)\ - K(t) ^-^LI) , (g.l) 



with the body frame components of the angular velocity being a;''(£"(T), £"(t)) = 

-ie™^ [i?^i?]"''(£"(T),£"(T)) . 

As already said, for every equation of state the mass density A(r, cTo) is a suitable function 
of no((?o)j K'^i'T^Bo) /n1{(Jo) and no{(fo)/ dct (^)- While the last term (absent only in the 
case of dust, because p = 0) is determined by Eq.(6.31), from Eq.(6.30) we get the following 
expression of the second term 

^^^ = ^2^ E rS(ao)r^,(a„)p,,(r)-p,,(r) = 

= ^27^ [ E i^o)) " (^) + 2 (a„) (a„) p^, (r) • p^, (r) + 

+ 2Er£(^o) E ^ni^o)Mr)-Mr) + 

i=l RJ^iLl 

+ E E rJ^(^o)rJ,(ao)p^,(r)-p^,(r)]. (8.2) 

By using Eqs. (7.10) and (7.31) this term can be expressed in the non-canonical basis 
, I TV I, ^, Ri, Hi, R2, ^n^ui.uj' Pn^ui,u2' '^^^ result is that the invariant mass density 

i) is independent from the Euler angles e^; 

ii) contains terms bilinear and linear in the body frame components 5^ of the spin. 

However, A(r, (To) is a complicated function of these three terms. The simplest expression 



is obtained in the case of dust, where A(r, Bo) = ylA* Wo(^o)]^ + Bo). Since the body 
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frame components of the angular velocity are defined [13,14,22] as 

Wrie^ = = F,,(5",g^p^)r + (8.3) 

it turns out that there is no linear relation between the spin and the angular velocity like in 
the non-relativistic rigid body (this property is true also for non-relativistic non-rigid bodies 
[13]). 

We can write the Hamilton equations for the orientational variables e^, and for the 
shape variables. Prom them we can deduce the equations of motion for the orientational 

2 — * 

variables a, S , 7 in Eqs.(7.10) (the other three | <S |, S^, (3 are Noether constants of motion). 
These three variables are not constant of motion for deformable bodies: they are coupled 
to the shape variables and describe how the dynamical body frame rotates when the body 
changes its shape. In particular a, being conjugate to the constant of motion | 5 |, is an 
ignorable variable {Ai, expressed in the canonical basis (7.10), does not depend on it). 

Three types of configuration for the motion of the fluid are interesting: 
i) Pure rotational motion - It is defined by constant shape configuration variables q'^ — 
0. With this condition the Hamilton equation q^^{q^,M.} become a system of algebraic 
equations for the shape momenta p^°^ — Piu\q=o- Even if we cannot find the explicit solution, 
its form is of the type p^°^ = Z^r ^ ^^(^) ?'')• The purely rotational Hamiltonian is At^'""*) — 

A4\ _ (o). Therefore for the dust invariant mass density we get 

p^—Pfji 



A(-*)(T,a„) = /[//n,(a„)]2 + $: A-(a„5M(r),5")5^5^ (8.4) 

y rs 

However Eqs.(8.1) show that the generalized Eulerian coordinates (i.e. the flux lines in 
adapted coordinates) do not perform a rigid motion: 

or ^ „ dpi' 

^{u){T)xtiT,ao)y. (8.5) 



^^Prom now on we shall use the following notations: i) q^^, p/j, will denote all the canonically 
conjugate shape variables \N\, ^, Ri, Hi, R2, II2, rn-tui,u^^ Pn^ui,u^'i ii) 1°"^ Pa will denote all the 
canonically conjugate shape variables p^j^ui including the first three pairs connected 

with the choice of the body frame axes. 
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ii) Pure vibrational motion - It is defined by the vanishing of the angular velocity, so that 
Eq.(8.3) determine the body frame components of the spin in terms of the shape variables: 
S^\uj=o — ^\o){q^tPij) 0- By putting this expression in A4 gives a purely vibrational 
Hamiltonian Al^"^^). 

iii) Small shape momenta - We can study configurations in which the shape momenta are 
very small. We can define the following two approximations (we use the dust to illustrate 
them) : 

3-) P/i ~ - The dust invariant mass density becomes 



A(t, Bo) ~ LnoiaoW + J2C'^i^o,q^)S''S'. (8.6) 

y rs 

— * 

If we make a Taylor expansion around the canonical center of mass Q = of the function 
[see also Eq. (9.2)] 

f-(T, ao) = det-^ j A(r, a,) (8.7) 

and if we define 

we can define the following new A- multipolar expansion of the invariant mass {V{t) is the 
volume of the fluid) 



M ~ V{r) Do{r) /[// n„(a„)]2 + J] C-(0, qi^) 6^5' + 

y rs 

+ E E Tr...rStX) I ^ d'cjo det — (8.9) 

n=l n-r„ JViO) \OOoJ 



where 

T 



The first term in this multipolar expansion is just the relativistic rotator defined in 
Eq.(5.12) of Ref. [27] and C^*(0, 5^) plays the role of the inverse of the tensor of inertia. 

b) Only Pa ~ - If we denote [/* the 3 shape variables |A^|, R2 and Vi their conjugate 
momenta ^, Hi, 112, the dust invariant mass density becomes 
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where if(3) is the non-relativistic vibrational (namely only function of the shape variables) 
Hamiltonian for the 3-body problem [13]. If we repeat the previous Taylor expansion (A- 
multipolar expansion), the first term will be 

V{t) D,{t) Lno{ao)f + Y.C''{^. U^)S'S' + ^ i?^'(0, g", U\ V,) + i7(3)([/\ V^) 

V rs r 

It corresponds to a generalized rotator interacting with a 3-body problem parametrized by 
the two dipoles i?2, originally along two of the chosen body frame axes, and by the angle 
between them, described by |A^|. 
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IX. DIXON'S MULTIPOLES. 

In this Section, following the general treatment given in Ref. [17], we shall give Dixon's 
multipoles [16,17] for the energy-momentum tensor of the perfect fluid in the rest-frame 

— * 

instant form on the Wigner hyper-planes in the gauge T t and Q^Q. 

On the Wigner hyper-planes the energy-momentum tensor (2.22) is rewritten in the form 



T^S(r, a) = I d'ao det 



(5^(a-E(r,a„))r^^(r,a<,), 



with [U^ and TZ are defined in Eqs. (3.16) and (3.10), respectively) 



nAB, 



U^{t, Bo) 



dp 
dh 



{r,Bo), 



1 (9S(r, ao) x 

7^(r,a„)V' dr 



f^^{T,ao) = A{T,ao) det 



-1 / di: 



daoj ' 



f^'-{T,ao)^K'{T,ao) def 



.1 / 9S 



dao ' 



f-(T,a<,) 



det"^ 













(9.2) 



Let us consider a world-line w'^{t) = a;^(r) + ^'^{p) (^'{t), where a;^(r) is the centroid, 
origin of the 3-coordinates on Wigner hyper-planes. Dixon's multipoles of the energy- 
momentum tensor with respect to this world-line are defined as 

^Mi...MnM<^(^(^)) = J £a (z^i(t, a) - w'^i(r)) ... {z^"{r, a) - w^-{r)) T^'^'izir, a)) = 



(9.3) 



The rest-frame instant form multipoles are 



g'"^"'^"^^(C(T)) = jd^'a {a^^ - C{t)) ... {a^- - C(r)) T^^r, a) = 



dV„ (S'-Ht, do) - r (r)) ... (E^H^, Bo) - r (r)) det 



d(Jo, 



T^^(r,a,). (9.4) 
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The monopoles are q'^^{C{'T)) with the mass monopole q^'^{C{'T)) — -M. and the vanishing 

— * 

(due to the rest frame condition) momentum monopole q'^^{C{T)) —V^^O. Then there are 
the stress tensor monopole q'^'^iCi^)) and the trace ^'^a (C (''")) • 

The dipoles are q^^^ {({'''))■ ^ we ask for the vanishing of the mass dipole 
Q^^'iCir)) = / d'ao [E^(T,a<,) -C(t)] A(T,a<,), 

we find 

5--(C(t))=0 ^C(t)=^~Q~0. (9.5) 

This means that the vanishing of the mass dipole identifies the world-hne of the internal 
M0ller center of energy, namely the centroid a;^(r), which, in the rest-frame instant form in 
the gauge Q ~ 0, is also both Tulczyjew and Pirani centroid, as shown in Ref. [17] in the 
case of particles. 

Therefore the multipoles with respect to the center of energy are 

= / dV„ (E^i(r, a„) - 7^^^T)) ... {^''{r, Bo) - n^'{r)) det f^^{T, a,) = 

= E ^S.(T)..r-(r)/dV, (rJ^(T,a„)-/i^,(r))...(r^„(r,a„)-/i^„(r))x 



mi,..,rnn 



X 



det lE^T^-nMl T^^(r,ao). (9.6) 



The last lines give the expression of Dixon's multipoles in terms of the relative variables, once 
the energy momentum tensor (9.2) is rewritten in terms of them. Then we could re-express 
the multipoles in terms of the orientational and shape variables. 

Then the momentum dipole takes the following form 

g--(T) = I d^Go (S'-(t, Go) - n^{r)) K%T, a^) = 

- / d'ao E rUr) (r^(T, a,) - hj^r)) K^{t, a,) ^ 

^e^;(^)p:M' (9-7) 



57 



and the angular momentum with respect to the internal center of energy, coinciding with 

— * 

the internal spin S, is 



As shown in Ref. [17] the quadrupoles q^^^^^^{T) allow to introduce two definitions of 
barycentric tensor of inertia: 



i) Dixon's one using the mass quadrupole, Idlllni'^) — ^^^^^ ?""^^(''') — g^^^^"^"^ (r) ; 
n) Thome's one, F^^^Ar) = 5''''' E„ ?""^a(t) - ^^^^^^^(t); 



both definitions give the standard tensor of inertia in the non-relativistic limit, since their 
difference is at the post-Newtonian order. 

In Ref. [17] there is the study of the non-relativistic limit of Dixon's multipoles by 
means of the Gartenhaus-Schwartz transformation, which are then compared with the non- 
relativistic multipoles defined in Appendix A of that paper. 





(9.8) 
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X. CONCLUSIONS. 



We have developed a new Lagrangian and Hamiltonian formulation of perfect fluids in 
terms of generalized Eulerian coordinates E(t, Bq) , which not only produces less complicated 
invariant masses, but also allows to define dynamical body frames, spin frames, orientational 
and shape variables by a natural extension of the techniques developed for N-body systems. 
On the contrary, it is too difficult to determine these quantities in the description based on 
the Lagrangian comoving coordinates. Since the Lagrangian contains the fluid density at 
the initial time, the choice of the two first axes of the dynamical body frame can be adapted 
to the initial form of the fluid. Then the dynamical body frame will evolve in time according 
to the changes in the form of the fluid according to the initial data and to the equation of 
state. 

We have also evaluated Dixon's multipoles in the rest-frame instant form. While they 
allow a description of the mean motion of the extended body, the fluid in this case, the 
orientational and shape variables give a complete information about the real motion with 
its changes of shape, moreover adapted to all the generic Nocthcr constants of motion. The 
invariant mass of the fluid, i.e. the Hamiltonian governing the real motion, can in turn be 
expressed with a A-multipolar expansion, as shown in Section VIII, which may be more 
useful than Dixon's multipoles when only a flnite number of shape variables is relevant (the 
others may be treated as perturbations) . 

We hope that this description of the fluid as an extended deformable relativistic body 
will help to treat with numerical simulations and/or approximations any type of system from 
droplet models of the proton, to heavy ion flreballs, plasmas and rotating stars. Regarding 
the simulation of stars we still need to make the coupling to tetrad gravity and then to go to 
a completely flxed 3-orthogonal gauge with the technique developed in Ref. [20] in absence 
of matter. This will give a new starting point for Hamiltonian numerical gravity and for 
the simulation of the properties of rotating stars. In particular in the case of incompressible 
fluids it will be interesting to try to recover the ellipsoidal equilibrium conflgurations in the 
non-relativistic limit. In the meantime in Ref. [22] it will be shown that the non-relativistic 
limit of the formalism of this paper allows to describe such conflguration after the addition 
of the Newton gravitational potential. 
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APPENDIX A: DEFINITIONS AND RELATIONS ON THE SPACE-LIKE 

HYPER-SURFACE 



The vierbeins z'^{T,a) = ^^l^'^f^ satisfy the Gauss- Codazzi-WeingaHen\niegYab\\\iy re- 



lation 

^4(r,g)--£j4(r.g) = o. (Ai) 

By construction the three 4- vectors Zf'{T, a) (f = 1, 2, 3) define a non-orthonormal basis 
for the tangent space to hyper-surface S(t); then it is possible to define the induced metric 
on the hyper-surface 

9fs(r, a) = z^(r, a)r]^^z^^{r, a), (A2) 

and the normal unit vector 

P(t, a) = 6^"^^^i,(t, a)z2p{T, a)z,,{T, a), (A3) 



where 



so that 



7(r,a) = -det(^,,(T,a)), (A4) 



r(r,a)/^(T,a) = l. (A5) 



Equally if we define 



g (r, a) = - det {g^B ir,a)), ( A6) 



we have 



UT,a)z>^{T,a)^ 



9[r, (y) 



We can define the inverse vierbeins z^{t, a) such that 



(A7) 



^iir, a)z't{T,a) = Si, 

z^JT,a)zj{T,a)^r]i:. (A8) 
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Then the inverse metric g^^{T, a), such that 



(A9) 



is defined by 



(AlO) 



We have also 



r7'^^ = 4(r,a)/^(r,a)4(r,a). 



(All) 



It is useful to consider the inverse 3- dimensional metric 7''*(t, a) such that 

r{T,a)grs{T,a)^5l. 



(A12) 



The following relation holds 



(A13) 



Prom this expression we can find 



If we define the lapse function 



N{T,a) 



and the shift function 



(A14) 



(A15) 



(A16) 



we can write the following expression for the inverse metric 



/^(r, a) = Y'ir, a) + — 1— iV^(r, a)N%T, a) 

M [r, a) 
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(A17) 



and we have 

^^(r, a) = 7V(t, a)F{T, a) + N\t, a). (A18) 

Moreover from the definition 

g^rir, a) = <(r, a)r/^.<(r, a), (A19) 

we have 

Qrrir, a) = N'ir, a) + N'{r, a)grf{T, a). (A20) 
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APPENDIX B: PERFECT FLUIDS ADMITTING A CLOSED FORM OF THE 

INVARIANT MASS. 



Let us rewrite Eq.(4.8) in the following form 



dX J 



(r, ao) X' - (r, a„) = A{t, a„) X\t, ao), 



B{T,ao) 



det (ft) 



A{t, ao) 



(Bl) 



For each equation of state p = p{n, So) with n = X the solution of this equation allows 
to get the explicit phase space form of the constraint 7i^(r, a) ^ in Eqs.(4.9). Referring 
to Section 5 of Ref. [6] for the determination of the equations of state p = p{h, So), in this 
Section we will show the few cases in which the solution for X{t, Bo) can be obtained in 
closed form. 

1) As shown in Ref. [6], for the dust we have p = and p{ti) — /in with /i — const. By 
using Eqs. (3.10) and (3.11), the action (3.15) becomes 





1=1 



(B3) 



Since we have p{X) - X ^ 0, Eq.(4.8) has the solution 



In Ref. [6] X = ^J^n was used as unknown. 
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and the second of Eqs.(4.9) becomes 

^±(T,a)=p^(T,a)P(T,a) + 



(B4) 



- y _ S(r, c7„))7[/ino(ao)]2 - 7-(r, S(r, a,)) ;^,(r, (T„) i^,(T, a^). (B5) 

On space-like hyper-planes and on Wigner hyper-planes the invariant mass and the in- 
ternal boost of Eqs.(4.27), (5.21), (5.23) become 



iC^ - j d^Go E(t, Bo) \J\pi-no{Bo)Y + K'^{t,Bo)- 



(B6) 



2) Let us now consider some cases of barotropic, p — p{p{n, §0)), and isentropic, 
P — P(p(^))) fluids. Let us remember that the dominant energy condition on the energy- 
momentum tensor requires \p\ < p. 

2a) p — kp{h) {k ^ —1), whose equation of state is p{n) — p,n^^^. With X = n we get 
^^{k + \)iiX^ and Eq.(Bl) becomes 

{{k + \)pfX'^''-'-\x'-B')=A. (B7) 

This equation can be solved in various cases, in particular for A; = 1 and k — \ 

A) k = 1, p = p, p = pn^i we have: 



A 



(B8) 



and then 

^±(r, a) = p^{t, a) l^{r, a) - f d'ao S\a - E(t, Bo)) 



2no{Bo) p 



^/^det(£) 



+ 



+ y^det( — )/i 



9S\ lY\^)Kr{r,Bo)Ks{T,Bo) , nl{Bo)p 



+ 



Vn.(a„) ^det^ (g) 



0, 



M 



d^Un 



2nJBn 



+ det 



9S 



, T, (Jn 



(B9) 
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3' -f^ 3 



B) Photon gas, k — ^, p — |p, p — jin^^^. Eq.(Bl) is the following cubic equation in 



If we define 

6 

3 r" ' \'ip 

to find the solution of the previous equation is equivalent to solve the cubic equation 



Z = Y-- Ub^+(—\ am , (Bll) 



+ Ci Z - Co = 0, (B12) 

where 

Using the Cardano solution we obtain: 

A 27 A\An) 27 \An ) 

Co , fT,\^ , [Co 



and then: 



2 +lf"^-°' ■ 



We can obtain the explicit expression of the constraint 7Yx(t, a) ~ and of the invariant 
mass M. using the solution (B15) in the eqs.(4.9) and (4.15). Since we have 

^4 ^2 I dp A 2 4 1 

= - p = ^/.xt = ^/.yi, (BI6) 

we obtain {Y = X^) 
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(BIT) 



2b) A variant is the equation of state p{n) = mh + ■^^{mhy' (7 7^ 1) with p = 
k (m n)'^ — {j — 1) {p — mn). It is called a polytropic fluid (7 = 1 + ^) by some authors and 
we can have k — k{so) in the non-isentropic case. Since we have 

dp 



dX 



m 



1 + ^— A;(mX) 
7-1 



7-1 



Eq.(Bl) becomes 



1 + 



7-1 



-1 



^AX\ 



A) For 7 = 2 we have the fourth order equation in X: 



m^(X^-B^) l + 2kmX 



^AX\ 



(B18) 



(B19) 



(B20) 



B) For 7 = 3 we have a third order equation mY — X'^: 



{Y - B") 



1 + 



"ikrn^ 



Y 



^AY. 



(B21) 



3) Standard polytropic perfect fluids have p = k p^{n) (7 = 1 + ^ 7^ 1) and 

p{h) = mh 1 — k {mhy~^ 



7-1 



Since we have 



9p 
dX 



m 



1 — k (mh) 



7-1 



7-1 



Eq.(Bl) becomes 



m 



(1 - km^-^ X^-^) {X^ - B^) = AX' 



For 7 = 3 (n = I) it is a fourth order equation m.Y — X"^: 



m 



(1 - A;m2y)^(y-S2) = Ay. 



(B22) 



(B23) 



(B24) 



(B25) 
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APPENDIX C: ON THE POISSON BRACKET 



In Section III we observed that the constraints of the Hamiltonian formulation are not 
exphcitly known as function of the canonical variables. Nevertheless it is possible to calculate 
their Poisson Bracket with a another functional on the phase space. To see this, we define 
the following short notations 



B 



R = 





p-X 


dp 






dX 


i( 


p-X 


dp 






dX 


X 


m 


-1 


2B 




■ 



(CI) 



With this notation the implicit definition of X, Eq.(4.9) can be rewritten in the form 

2 



dp 

ax 



(C2) 



and the second constraint of Eqs.(4.10) is 

W±(r,a)=p^(T,a)P(T,a)-^(T,a), 

where 

T{T,a) = J d^aoUoiao) - E(r, ao))Q. 
From Eq.(C2), using the rules on the Poisson Bracket, we can see that 

dp 



(C3) 



(C4) 



+ 2 



dX 



d^p B^ d^p df^B^ 
"dX^X^ ^ dJO ^ dXX^ 



{.,X}. 



(C5) 
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Equally we can see that 

{., ^(r, a)} = -l d'ao —5\a - S(r, a,)) Q {., S^(r, + 



+ 



X 



j (faono{Bo)5\a-t{T,ao)) x 



B'^V dX X dx 



X 



s I dx^ X dx^ x^ dx 



{.,x}. 



Using the results (C5) we can write 



1 ^] I B d'p 1 dp 

B\ dX^j X ^2 9^ 



2S \dX 



-1 



{.,A} + 2i? 



axi x^ 



(C6) 



(C7) 



If we substitute this expression in Eq. (C6) and using the short notation (CI) we obtain 

d 



+ 



+ 



~UB}-R{.,A} 



(C8) 



Being A and B known functions of the canonical variables, the previous equation permits 

to calculate the Poisson bracket with the constraint 'Hj_{t,(t) although the X is unknown 
explicitly. The rule (C8) can be used, for example, to verify the algebra (4.13). 



On the hyper-planes, the previous observations are again valid; in this case we have 



B = 



(C9) 



nlia.) ' det(g)- 

From Eq.(4.27) the invariant mass and the canonical generator of the internal boost are 
defined in term of the density 
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A(t, ao) 



det 



'as 

dan 



dp 



1 dp 



such that 



On them the previous rule becomes 



{., A(t, (To)} = no(ao) 



-^{■,B}-R{.,A} 



(CIO) 



(Cll) 



(C12) 



This rule allows us to calculate the Poisson brackets of a functional on the phase space 
with the constraint also in the case of hyper-planes or Wigner hyper-planes. This is useful 
for example for getting the equations of motion from the Hamilton-Dirac equations. 
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APPENDIX D: GARTHENAUS-SCHWARTZ TRANSFORMATIONS 



Using the notations of Section VI, let G — V ■ Qhe the canonical generator of a trans- 
formation, called Garthenaus- Schwartz canonical transformations [29]. If F is an arbitrary 
functions on the phase space, we have that its infinitesimal transformation is 



SF ^Sa- {F,V ■ Q}. 

For finite values of the parameter a the transformation can be written as 

F{a) = F+ r da {F{a),V{a) ■ Q{a)}. 
Jo 

We are interested to the singular limit a — > oo and we use the notation 

F' = lim F{a). 



[Dl) 



(D2) 



(D3) 



Deriving both sides of (D2) we can see that to realize the canonical transformation is 
equivalent to solve the differential equation 



It is trivial to verify that 



= {F(a),V{a).Q(a)}, 
F(0) = F. 

Via) = e-'^V ^ 



(D4) 



(D5) 



The usefulness of the singular hmits is shown by the following observation [29] . Let F 
be a function on the phase space such that 



{P,F} = {P{a),F{a)}^0. 



(D6) 



For this function the singular limit 



lim F(a) = F', 



(D7) 
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exists and is well defined. Moreover, if we define 

G^{Q,F}, (D8) 

the Jacoby identity implies that 

{V,G}^0, (D9) 

and the limit 

lim G(a) = G', (DIO) 

a— >oo 

exists and is well defined. In conclusion 

{r, F'} = lim {V, F{a)} = lim e+°{:P(o;), F{a)} = 0, (Dll) 
because {V{a), F{a)} = 0, and 

{Q,F'}= lim{Q,F(a)}= lim e-"{Q(a), = 0, (D12) 

because the singular hmit of {Q{a), F{a)} is the well defined quantity G'. 

This observation can be applied to the relative variables W{t, ao), P^{t, ag). These vari- 
ables are by construction such that 

{V, r (r, ao)} = {P, p'{r, a,)} = 0, (D13) 

and then their singular limits 

Jim r(a I T,ao) = 3fJ"'(r,ao), 

lim p'{a I r, ffo) = p"{t, ffo), (D14) 

a— +0O 

exist and are well defined. In particular we get 

{V,^'^{T,ao)} ^ {V,p'^{r,ao)} = 0, 

{Q,^'^(T,ao)} = {Q, p"-(r,a,)} = 0, 

{5ft"- (r, ao), p'^ir, a'J} = 5" 5'{a, - a',). (D15) 
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In conclusion the coordinates of the table 









Q 






V 







are canonical coordinates. Moreover we can observe that 



^"{T.do) = 3?''(r,ao) + / da {W\a \ r,ao),V{a) ■ Q{a)} 

Jo 



Jo 

/'■(r, Bo) = p''(t, Bo) + da {p''{a \ r, Bo),V{a) ■ Q{a)} = 

Jo 



P'ir, Bo) + / da e-"{p^'(« I ^, ^o), Q\oc)} 
Jo 



^pr^r,Bo)+Vni\ 



(D16) 



(D17) 



Due to the previous considerations, the {W{a \ T,Bo),Q.^{a)} and the {p^{a \ 
'T,Bo),Q^{a)} are well defined constants in the singular limit. Then the integral I^,Ip 
in Eq.(D17) are well defined for the presence of factor. Then if we use exphcitly the 

— * 

condition we get 



3?"^(t, Bo) ^ r (r, Bo) p''{t, Bo) ^ p'{t, Bo). 
This is the case of the gauge fixing (6.27). 



(D18) 
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APPENDIX E: SOME SOLUTIONS FOR THE KERNEL F 

In this Appendix we construct some kernels F that satisfy the conditions (6. 5), (6. 8) and 
(6.10) or Eq. (6.18). 

The first solution is based on the possibihty to read the kernel as distributions. For 

example, we can define in this case 

rK{ffo-^'o)-^lJ'{^o-K), (El) 

and the second of Eqs. (6.5) is satisfied being reduced to 

Vi 1 = 0. (E2) 

Let us define the usual symmetric Green function c{ao — a'^) 

VlA^o-o':)^5\ao-a':), (E3) 
and let us make the following ansatz on F^ 

rE(a., a;) =c(a„ -<)+/«). (E4) 
Then Eq. (6.10) is satisfied 

Fs(a„, a;) = c(a<, - <) = 5\ao - a',). (E5) 

The function / is determined imposing Eq.(6.5) and we get 

fi^'o) = - j^ J c?Voino((?oi)c(aoi - a^). (E6) 
Automatically also Eq.(6.8) is satisfied. In conclusion 

rs(^o, ff'o) = c{ao ~^'o) ^ Jf J no{aoi)c{aoi - Bo), (E7) 
is a distribution-hke solution for the kernels F. 

Another class of possible solutions is obtained if we use the representation F^((To), r^(ao) 
given by Eq. (6.17). Then we can consider a second base of orthonormal functions ^„(ao) 
such that 
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Jd^aono{ao)'ifo{^o)y^O. (E8) 

The base no{^o),'^nyto{^o) is a complete, non orthonormal base of functions. Using the 
Gram-Schmidt procedure [34] we can construct the orthonormal base "^'^{ao) such that in 
particular 

Ki^o) = (E9) 

with the normalization constant 

J d^aonliao). (ElO) 

The other elements of ^^^o((?o) ai'c given by the recurrence formula of the Gram- 
Schmidt's algorithm. With these definitions the first of the conditions (6.18) is satisfied 
if we choose 



rl{ao)^K{^o) if n^o. (Ell) 
The fourth of Eqs.(6.18) is satisfied if the F^'s have the following form 

rf (a'o) = 0, 

T^{ao) = K{ao)-CnK{^o) if n + o. (E12) 



We use the second of Eqs.(6.18) for fixing the values of the coefficients c„ 

R 



c„, = -— j d^Oo <(ao) n + o. (E13) 



With this choice also the third of Eqs.(6.18) is satisfied. In fact we can calculate explicitly 
the sum in the left-hand side of this conditions using the completeness of the basis ^'^(ao) 

= ^\o^o - B20) - K{a2o)Kia,o) + E %*;(^2o)*;(^io) = 

= 5\B^o - + <(a2o) E Cn<(^lo), (E14) 
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with Co — —1 in the last hne . Finally we can observe that, in the sum, the c„ are the 
components on the base ^^(^o) of the {—R/M) constant function and then we have in 
accord with Eqs.(6.18) 

K{^2o) E C^K{B.o) = -%{^2o)f^ = (E15) 

In conclusion 

rs(<?lo, ?2o) = E *;(<?lo)^n(<?2o), 

TK{a,o. a2o) = E ^n{a,o)K{^2o) ' E % ^n(^lo)- (E16) 
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